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: 1 Introduction 
cn 

Q ■ Consider a dynamical system 

^ ' together with its stochastic perturbations 

: dXf = b{X^^')dt + ea{X^'')dWt, Xo"'" = x G M'^. (2) 

H : 

■ Here e > is a small parameter, Wt is a Wiener process in R'', and the coefficients a and b 
are assumed to be Lipschitz continuous. The diffusion matrix a{x) = (ajj(x)) = a{x)a*{x) 
is assumed to be uniformly positive definite. 

Together with ([2]), we can consider the corresponding Cauchy problem 

u%0,x) = g{x), X G (4) 
where 5^ is a bounded continuous function. 
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Suppose for a moment that the vector field b has just one asymptotically stable equi- 
librium point O such that all the points get attracted to O and {b{x),x — 0)< —c\x — 0\ 
for some positive constant c and all sufficiently large Then it is easy to check that 

lim F(Xf'' eU) = l 

(£,t)^{0,Oo) 

for any neighborhood U of the equilibrium O. Taking into account that the solution of 
([3])- (HI) can be written in the form u^(t, x) = Egi^X^'"^) and the continuity of g, we conclude 
that 

lim -u^ft, x) = qiO). 

A similar result holds in the case of a unique compact global attractor if the system ([1]) has 
a unique normalized invariant measure on the attractor. This is the case, for example, if 
system ([T]) in has a unique limit cycle attracting all the trajectories except the unstable 
equilibrium inside the cycle. 

The situation becomes more complicated if the dynamical system has more than one 
asymptotically stable attractor. Assume, for brevity, that all the attractors are equi- 
libriums Oi,...,0„. Let Di be the basin of Oj, 1 < i < n, and assume that the set 

\ {Di U ... U Dn) belongs to a finite union of surfaces of dimension d — 1. The long time 
behavior of X^'^ and u^{t,x) is now determined by the transitions between the attractors 
Oi,...,0„. These transitions are described by the large deviation theory for stochastic 
perturbations of dynamical systems developed in the late 1960-s (see [7j and references 
there). In particular, the weak limit ^ of the invariant measure //^ of the family of pro- 
cesses (|2]) was found. In the generic case, the limiting measure /i is concentrated on one 
of the attractors, which will be denoted by O*. Then 

lim lim u^{t,x) = g{0*). 

However, in the case of many attractors, the limiting behavior of Xf'*^ and u'^{t,x) as 
e I and t ^ oo depends on the way in which {e,t) approaches (0,cxd). Roughly 
speaking, under natural additional assumptions, there exist a finite number of regions 
in the neighborhood of (0, oo) such that the limiting distribution of ''^ and the limit 
of u'^{t,x) exist if {e,t) approaches (0, oo) while staying inside one region. For different 
regions, these limits are, in general, different. 

The corresponding theory of metastability (of sublimiting distributions) was developed 
in [3] (see also [3], [7], [S]). The notion of a hierarchy of cycles, which is discussed 
below, was introduced there. Let S'o,t(v^) be the action functional for the family X^'^ in 
C([0,T],R'^) asejO ([?]): 

SoA¥>) = ^ r E ('''(VtM - Hvt)m - bj{^t))dt, T > 0, ^ G C([0,T],M'^), 

"^0 i,j=l 
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where a^^ are the elements of the inverse matrix, that is a^^ = (a and S'o,t(v^) = +00 
if ip is not absolutely continuous. The quasi-potential is defined as 



Note that while the term "quasi-potential" is normally applied to the function V of the 
variable y with x being a fixed equilibrium point, we use the same term for the function 
of two variables. The hierarchy of cycles is determined by the numbers 



The equilibriums Oi, ...,0„ are the cycles of rank zero. In the generic case, for each 
there exists a unique "next" equilibrium Oi = M{Oi) defined by Vu = min^^j Vik- For 
all sufficiently small 6 > 0, with probability close to one as e | 0, the process that 
starts in a (5-neighborhood of Oi will enter a 5-neighborhood of N'{Oi) before visiting the 
basins of any of the equilibriums other than Oi and Af{Oi). The time before the process 
enters the neighborhood of Oi = MiOi) is logarithmically equivalent to exp(Vi//5^). If the 
sequence O^, 7V(0i), M^{Oi) = M{U{Oi)), ...,W{Oi), ... is periodic, that is A/'"(Oi) = 
Oi for some n, then a cycle of rank one appears. It contains the cycles or rank zero 
Oi,U{Oi), ...,Ar"-i(Oi). If A/'"(Oi) ^ Oi for any n > 1, we say that Oi forms a cycle of 
rank one. The entire set of equilibriums is decomposed into cycles of rank one, which will 
be denoted by C\, C^_^. Note that some of the cycles of rank one may consist of one 
cycle of rank zero. 

Next, the transitions between cycles of rank one can be considered. Namely, in the 
generic case, for each cycle Cj there is a different cycle N'{C}) of rank one determined 
by Vij, 1 < i,j < n, with the following property: if the process starts at one of the 
equilibrium points in C}, then, with probability close to one as e | 0, it will enter a 
^-neighborhood of one of the equilibrium points inside the cycle M{C}) before visiting 
basins of any of the equilibriums outside C} and MiC}). This leads to the decomposition 
of the set of cycles or rank one into cycles of rank two. 

This procedure can be continued inductively until we arrive at a single cycle of finite 
rank R which contains all the equilibrium points. The cycles of rank r < R will be 
denoted by C[, ...,C;;^^. 

Let T^(A) = exp(A/£:^). (The results stated in the paper also hold for T^(A) x 
exp(A/£:^), that is if InT^(A) ~ A/e^ as e J. 0.) In the generic case, there is a finite 
set A C (0, 00) such that for each x E Di U ... U Dn and each A G (0, 00) \ A, one equilib- 
rium Om(x,x) is defined such that the measures fi^(T) = P(X^f^^^ G F) converge weakly to 
the (5-measure concentrated at Om(x,\)- The state Om{x,x) is called the metastable state 
for the initial point x and the time scale T^(A). 

In this paper, instead of the linear problem ©-(jl]), we will consider the Cauchy 
problem for the quasi-linear equation with a small parameter 



V{x, y) = inf {^o,t(<^) : V ^ C{% T], W"), ^{0) = x, <^(T) = y}, x, y G 



Vij = V{0i,0j), l<i,3<n. 



du'^{t, x) 
di 



.2 



d'^u'^{t, x) 
dxidxj 



L'u' ■= — 




+ b{x) ■V^u'{t,x), xgM"^, t>0, (5) 
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u'{0,x) = g{x), X e R'^. (6) 

We assume that the coefficients of this equation are Lipschitz continuous and bounded; 
the matrix {aij{x, u)) is assumed to be uniformly positive definite. Under these conditions, 
problem ([5])- ([6]) has a unique solution for any continuous bounded g{x) (see, for instance, 

0)- 

For if: > and x G M^, we can define X*'^'*^, s G [0,t], as the process which starts at x 
and solves 

dX'f'^ = b{Xl''''')ds + ea{X'f'', u%t - s, X'f''))dWs, s < t, (7) 

where the entries aij, 1 < i, j < d, of the matrix a{x, u) are Lipschitz continuous and such 
that aa* = a. The process X^'^'^ can be viewed as a nonlinear stochastic perturbation 
of ([I]). As in the linear case, we have the following relation between m*^ and the process 

u%t,x) = Eg{Xt''n- (8) 

Conversely, ©-([HI) can be viewed as a system of equations with unknown function 
and process X*'^'^. Under the above assumptions on the coefficients and the function g, 
the solution of this system is unique and the function -u^ is the solution of the problem 
([5]) -([6]). This approach of examining systems of the type ([71)-([H]) was employed in [4J in 
order to study linear and non-linear degenerate equations. The ffist initial-boundary value 
problem for quasi-linear parabolic equation with a small diffusion and the exit problem 
for the corresponding processes were studied in [6]. The results of the latter paper will 
be used here. 

In the nonlinear case, the role of Vij, 1 < i, j < n, will be played by the functions 
V,,{c) = inf{5oV(<^) • e C{[0,T],R^),^{0) = 0„^{T) = O,}, c G M, 

where 5*0 ^(v?) is defined using the function a(x, c) instead of a{x) used in the linear 

case. If A is sufficiently large, then the distribution of X'^^j^^^'^'^ , even in a generic case, 
converges not necessarily to a 5-measure concentrated at an equilibrium point, but to 
a distribution on the set of equilibrium points. This happens, for example, in the case 
of two equilibrium points if Vi2(c) = V2i(c) for some value of c. Therefore, in the case 
of nonlinear perturbations, the notion of a metastable state should be replaced by the 
notion of a metastable distribution. Note that metastable distributions arise also when 
perturbations of nearly-Hamiltonian systems are considered (see [1], [2]), but because of 
different reasons. 

An additional difficulty in analyzing the asymptotics of the process and the solution 
to the parabolic equation is due to the fact that the entire hierarchy of cycles may change 
with time in the non-linear case. An example of this phenomenon is provided in Section [61 

In Section[2lwe introduce some of the definitions and discuss the notion of the hierarchy 
of cycles in more detail. In Sections [31 and [H we consider a system with two equilibriums 
and a system with three equlibriums on the real line in the case when the hierarchy of 
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cycles is preserved. In Section [5] we formulate a general result for the case when the 
hierarchy of cycles is preserved. In Section [6] we study the asymptotics of the solution 
to the parabolic equation for a system in which a bifurcation in the hierarchy of cycles 
occurs. 

2 Preliminaries and notations 

Let a{x) be a symmetric dx d matrix whose elements aij{x) are Lipschitz continuous and 
satisfy < Ylt j=i^iji^)^i^jy ^ ^ ^ ^ ^"'^ Let a*-' be the elements of the inverse 
matrix, that is a*-' = {a~^)ij, and o" be a square matrix such that a = aa*. We choose a 
in such a way that aij are also Lipschitz continuous. 

Let SQrp be the normalized action functional for the family of processes X^''^ satisfying 

dX^'' = b{X^'')dt + ea{X^'')dWt, X^'' = x, (9) 

where 6 is a Lipschitz continuous vector field on M'^. Thus 

for absolutely continuous ip defined on [0,T], ip^ = x, and S'^j.(v9) = cxd if is not 
absolutely continuous or if <^o 7^ ^ (see [7]). Let ^"(x, y) be the quasi-potential for the 
family Xf'^ in R'^, that is 

V^"(x,i/) = inf{5oX</^) • ^ e C([0,T],M'^),y.(0) = x,^(r) = y}, x,y e M'^. (10) 

Let V^'j = V^{Oi, Oj). For a given function a, we define inductively the following objects 
(see [3], [Zj for a detailed exposition). 

(a) The hierarchy of cycles CJ', C^^, r < R. 

(b) The notion of the "next" equilibrium v{Cl) and the "next" cycle N'{Cl) of the 
same rank for a cycle C[ of rank less than R. 

(c) The transition rates V^tq., 1 < < 1 < j < Oj ^ C[, from a cycle to 
equilibriums outside this cycle. 

For r = 0, we define = {OA, V^o n — Assume that the cycles of rank r 
and the transition rates from those cycles to equilibrium points have been defined. We 
define Oj to be the next equilibrium after C[ if mm.j,o^^ci Vc'r,Oj is achieved at j. 

Assumption A. The minimum minj:o^^c[ Oj is achieved for a single value of j. 

We will write Oj = v{C'i) to express that Oj is the next equilibrium after C[. We say 
that the cycle C'l of rank r is the next after C[ if CI contains v{Cl). We will express 
this relation by writing C[' = N'{Cl). Starting from a cycle C[ of rank r, we can form 
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the sequence Cl^AfiCD^Af^iCl), ... by using the operation "next". If this sequence is 
periodic, that is C[ = N"'{Cl) for some n, then the cycles C[, A/'"~^(C[) form a cycle 
of rank r + 1. If C[ 7^ N''^{Cl) for any n > 1, then CI is said to form a cycle of rank r + 1. 
This way, the collection of all the cycles of rank r is decomposed in a non-intersecting 
union of cycles of rank r + 1. 

If C{', CJ form a cycle of rank r + 1, which will be denoted by F, we define 

^r^ = max VSr^^.^cr^) + mm - ^cl,.(cs,)), ^ T. (11) 

We can continue this procedure until we arrive at a single cycle of rank R. 

If r is a cycle, we define Dy = Ui.OievDi. As follows from [3] , [Z] , if the process starts 
in Dy, where F is a cycle of rank r < R, then with probability which tends to one as e J, 
it will leave Dp and enter a small neighborhood of z^(r) in time T{e) x exp(l^^^p^/£^). 

Note that the quasi-potential can be defined by ( ITOi) even if a has some discontinuities. 
We shall be particularly interested in the structure of the hierarchy of cycles and the 
exponential transition times for functions a which are of the form a = a{x, f{x)), where 
/ is constant on each Di. The reason for that is that the solution of ([5])- ([6]) is nearly 
constant inside each of the domains Df = {x E Di : dist{x, dDi) > 6, \x\ < 1/6}, 6 > 0, 
1 < i < ra, as follows from the following lemma. 

Lemma 2.1. For every positive Aq and 6 there is positive eo such that 

\u'iT%X),x)-u%T'iX),0,)\<S (12) 
whenever x G D^, e < and A > Aq. 

Proof. A similar lemma was proved in [6j, so we only recall the main steps here. The solu- 
tion at time r^(A) can be viewed as the solution of the same equation at time T'^{X) with 
initial data g replaced by the solution of ([S])-([SD at time T^(A) — T'^(A) (semigroup prop- 
erty). We use a priori estimates to show that the solution at time T'^(A) — T^(A) is almost 
constant in a small neighborhood of Oj. Instead of ([8]), we can use the representation 

u%T%\),x) = Eu%T%\) - r^(A),X™'^'^). 

Finally, we note that X^^^j!^' ' is in a small neighborhood of Oi with high probability if 
A is sufficiently small and x E Df. □ 

3 The case of two equilibrium points 

In this section we assume that there are two asymptotically stable equilibrium points 
Oi,02 G W^. Let Di C be the set of points in which are attracted to Oi and 
D2 C the set of points attracted to O2. Note that in the case of two equilibrium 
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points, the hierarchy of cycles is always the same: Oi and O2 are cycles of rank zero, and 
there is one cycle of rank one which contains both Oi and O2. 

Let grain = '^'^^x&L'i 9{.^) a^d (7max = sup^-gj^d ^'(a;) . Define the functions M12, M21 : 

[5'min, fi-max] ^ ^ via 

M12{C) = M21(C) = r«t,Or 

These functions are shown on Figure 1. It is not difficult to check that the constant c in 
the second argument of a can be replaced by any function which equal to c on Di in the 
definition of Mu and equal to c on D2 in the definition of M21. 




Figure 1: The case of two equilibrium points 



Without loss of generality we may assume that g{Oi) < g{02). Let Ai = Mi2{g{0i)) 
and A2 = M2i{g{02))- In order to formulate the results on the asymptotics of u^(T^(A), x), 
we need the functions c^(A) and c^(A) defined as follows: 

= / siOi), < A < Ai, , . 

"^^^^ \ min{(7(02),min{c:cG [(7(Oi),(7(02)],Mi2(c) = A}}, A > Ai, ^'^^ 

r^(\) = S 9(0^^^ 0<A<A2, . . 

''^'^> \ me.x{g{Oi),max{c:ce[g{0,),g{02)],M2i{c) = X}}, A > A2. 

Let A3 = inf{A : c^(A) > c^(A)}. Assume that at least one of the functions and 
is continuous at A3. Let c* = c^(A3) if is continuous at A3 and c* = c^(A3) otherwise. 
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Let c^(A) = min(c^(A), c*) and c^(A) = max(c^(A), c*). On Figure 2, the graphs of and 
are denoted by the thick and the dotted hnes, respectively. 
The asymptotics of u^{T^{X),x) is described by the following theorem. 

Theorem 3.1. Let the above assumptions be satisfied. Suppose that the function {X) is 
continuous at a point A G (0, oo). Then for every 6 > the following limit 

limu%T%X),x) = c\X) 

elO 

is uniform inx E D{. Suppose that the function c^(A) is continuous at a point X E (0, oo). 
Then for every 6 > the following limit 



\imu'{T'iX),x) = c\X) 

e[0 



is uniform in x E -Dg- 



Proof. The proof essentially relies on the results of [6], where a similar statement was 
proved for the case of the initial-boundary value problem with one equilibrium point 
inside the domain. We only sketch one of the main steps since it illustrates the ideas 
involved in the entire proof. Namely, we will show that if < A < A3 and c^,c^ are 
continuous at A, then 

limsup sup ■u^(T^(A),x) < c^(A) 

and 



liminf inf u\T'{X),x) > c^(A). 
Due to Lemma [2. H it is sufficient to show that 



limsupM"(T"(A),Oi) < c^(A), liminf u"(T"(A), O2) > c2(A). (15) 

Note that there is a positive Vq such that for every < 5 < there is £0 > such 
that 

\u\T'{X),x)-g{0,)\<5 (16) 

whenever x E Df , Q < e < Sq and 5 < X < vq. Indeed, f|T6l) follows from the representation 
([8]) of the solution in terms of the process ^'^'>'^'^ and the fact that belongs 
to a small neighborhood of with high probability if Vq is sufficiently small, as follows 
from the basic properties of the action functional (see |7|). 

If (1151) fails for a certain value of A, then due to continuity of the functions M^(t, Oi) in 
t, it follows from f|T6l) that for an arbitrarily small 5' > there are sequences i and 
A„ E [5\ A] such that 

u'{t,Oi) <c\X) + 5\ u%t,02)>c^{X)-5', T'-{5') <t<T'-{Xn), 
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and either 

u'"iT'"{Xn),0,) = c\X) + 5' 

or 

M^"(T^"(A„),02) = c'(A)-y. 

Assume, for example, that the former is the case. Take 6" G (0, 6') which will be specified 
later. Due to continuity of u'^"{t, Oi) in t, we can find a sequence /i„ G [6', A„) such that 

and 

u'"it,0,) e [c\\) + 5",c\\) + 5'] for t e [T^"(/i„),T^"(A„)]. (17) 

We can express u^"(T^"(A„), Oi) in terms of the process xj ^nd the solution at 

the earlier time T'^"{fin) as follows 

«-(T-(A„),Oi) = En- (T-(^„),xj::(t^^^--^^^ • (18) 

Since is continuous at A, there are arbitrarily small 6' > such that Mi2(c^(A) + 6') > 
Mi2(c^{X)) = A. Since A„ < A, a process starting at Oi and satisfying ([9]) with 

aa*{x) = a(x,u-(T-(A„),Oi)) = a{x,c\X) + 6')) 

will be in an arbitrarily small neighborhood of Oi at time T'^"(A„) — T'^"{fin) with proba- 
bility which tends to one when e„ | 0. This remains true if the constant m— (T^"(A„), Oi) 
is replaced by a function which is sufficiently close to this constant in Df, where 6 is 
sufficiently small. Therefore, due to f|T7j) and Lemma 12.11 we can choose 6" sufficiently 
close to S' so that -^T=„(Ar)5^n'(/t„) ^^^^ ^ small neighborhood of Oi with probability 
which tends to one when e„ | 0. With S' and 5" thus fixed, we let | in f[T5]) . The left 
hand side is equal to c^(A) + 5', while the right hand side tends to c^(A) + 6". This leads 
to a contradiction which proves that (fT5|) holds. □ 

Remark. If A > A3, then c^(A) = c^(A) = c*. It is possible to show that the limit 

limM"(T^(A),x) = c* 

is uniform in (x. A) G Bi/s x [A, 00) for each A > A3. Therefore, for each S > and A > A3 
there is £0 > such that 

\u%t,x)-c*\ < 5 
whenever e G {0,eo), x G Bi/s and t > T^(A). 

Let xf^^^'^''", s G [0,T^(A)], be the process defined in ([7]). As follows from the 
large deviation theory (see [7j), the distribution of the random variable X^^^^^j'^'"^ will be 
concentrated near the points Oi and O2. From Theorem 13.11 and the representation ([8]) 
for the solution, we obtain the following theorem. 
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Theorem 3.2. Suppose that g{Oi) < g{02)- If the function {X) is continuous at a point 
A G (0,00) and x G Di, then the distribution of the random variable X'^^j^^}^'^'^ converges 
to the measure ^\ = 016^^ + 02^x2 > where the coefficients Oi and 02 can be found from the 
equations c^(A) = aig{Oi) + 02(7(02), Oi + 02 = 1. 

If the function (P {X) is continuous at a point X G (0, 00) and x G D2, then the distribu- 
tion of the random variable converges to the measure /ig = o-i^xi + (12^x2! where 
the coefficients ai and 02 can be found from the equations c^(A) = aig{Oi) + 02(7(02), 
ai + 02 = 1. 

// A G (A3, 00) and x E D, then the distribution of the random variable xJj^^^^^'^'^ 
converges to the measure fi* = aifJ^i +(12^x2; where the coefficients ai and 02 can be found 
from the equations c* = aig{Oi) + 02(7(02), ai + 02 = 1. 



4 Three equilibrium points without changes in the 
hierarchy of cycles 

In this section we assume that there are three asymptotically stable equilibrium points 
Oi, 02,03 such that g{Oi) < g{02) < giO^). For ci, 02,03 G [fi'min, 5'max], let 

fci,C2,C3ix) = CiXDiix) + C2XD2ix) + CsXDsix), X G R'^ . (19) 

Recall the definition of the hierarchy of cycles from Section El We will assume that, for 
each choice of constants q G [(7min, (7max] in the function a = a(x, Z^, 02,03(2;)), Assump- 
tion A holds and Oi and O2 form a cycle of rank one. Consequently Oi, O2 and O3 form 
a cycle of rank two for each choice of the constants. Define 

Mi2{c) = V^% M2,{c) = V^% 

Let Ai = Mi2(^(0i)) and A2 = M2i(^(02)). Define functions and by ([13]) and CT. 
respectively. Let A3 = inf{A : c^(A) > c^(A)}. Assume that at least one of the functions 
and is continuous at A3. Let c* = c^(A3) if is continuous at A3 and c* = c^(A3) 
otherwise. Let c^(A) = min(c^(A), c*) and c^(A) = max(c^(A), c*), A < A3. Let A4 = 
Mr3(c*) and A5 = MsrigiOs)). 

Let us assume that Ai < A2 < A3 < A4 < A5 (see Figure 2). For A < A3, the behavior of 
the solution in Di and D2 is still governed by Theorem 13. 1[ For each A > A3, the value of 
M^(T^(A), x) will be nearly constant on Df U D2, and we can treat the cycle F = {Oi, O2} 
in the same way a single equilibrium was treated in Section O Namely, let 

""^^^^ = { mm{(7(03),min{c : c G [c* , giO^)], Mrsic) = A}}, A >~A4, 
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X 



I12 



g(0,) 



t 



9(02) 



g(03) 



Figure 2: A case of three equilibrium points without changes in the hierarchy of cycles 



Define Ae = inf{A > A3 : c^(A) > c^(A)}. Assume that A5 < Ae and that at least 
one of the functions (F and (? is continuous at \. Let c** = if is continuous 

at Ag and c** = e'[\^ otherwise. Define c^(A) = c^(A) = min(c'"(A), c**), A > A3, and 
c3(A) =max(c='(A),c**), A > 0. 

Having thus defined the functions c*(A), i = 1,2,3, for all A > 0, we can now state 
that for each A > such that c* is continuous at A and every 6 > 0, the limit 



is uniform in x G -Df . 

On Figure 2, the limits \imeiou^{T^{X), x), as functions of A, for x G D(, D2 and 
are depicted by thick, dotted and dashed lines, respectively. 

5 A general result for the case when the hierarchy of 
cycles does not change 

In this section we will suppose that, in addition to Assumption A, the hierarchy of cycles 
and the equilibrium points z^(r) for each cycle F of rank less than R do not depend on 




limu"(T"(A),x) = c\X) 
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the choice of constants q G [(^min, S'max] in the function a = a{x, Yl'i=i CiXA(^))- 

We will say that a cycle F is active for a given value of A > if ^r^i,(r) < ^- We 
will say that it is engaged if ^,^(r) ~ ^ passive if V^^^j.-^ > X. We will say that a 
cycle To is connected to a cycle F by a chain if there is a sequence of cycles Fi, ...,Tk 
and equilibriums Oi G Ti,...,Ok G F^, Ok+i G F such that F, are engaged or active and 
Oi+i = z^(Fj) for < i < A;. The collection of all the cycles that do not belong to F and 
are connected to F by a chain will be called the cluster connected to F. For each cycle F 
of less than maximal rank and c G [gmm, dmax], we define 

Mr{c) = V^^^'[^^ 

and, for A > and C2 > ci, define C(ci,C2,A,F) = min(c2,inf(c > ci : Mr(c) > A)). 
Similarly, if C2 < ci, define C(ci, C2, A, F) = max(c2, sup(c < ci : Mr(c) > A)). 

In Figure 3 we have an example of a hierarchy of cycles with the thick arrows between 
the actively connected cycles and the corresponding equilibrium points. The dashed 
arrows are used for the engaged cycles and the dotted arrows for the passively connected 
cycles. 




Figure 3: The hierarchy of cycles 



In order to describe the asymptotics of m^(T^(A), x), we will define a finite number of 
"special" points = Aq < Ai < ... < A^ = oo. We claim that there are functions c*(A), 
1 < i < n, which are continuous on each of the intervals (A^, Afe+i), < k < m — 1, have 
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one-sided limits as A approaches the end points of the intervals, and are such that the 
limits 

\iuiu'{T'{X),x) ^c%X) 

are uniform in x E Df for each S > 0, X E \ A. Moreover, neither of the cycles 

changes its type (between passive, engaged and active) for A G {Xk,Xk+i) and a{x) = 
lim^io ■u^(T^(A), x)). We will use induction on k in order to define the functions c*(A) 
and describe for each cycle whether it is passive, engaged or active for A G (A^, A;^.+i) with 
a{x) = \ims^Qa{x,u^{T^{X),x)). In the process, we will make several assumptions about 
the functions Mr. 

Assuming that we have defined c*(A), let 

Ar = inf{A > : c* (A') docs not depend on i for A' > A and Oi G F}. 

From the inductive construction of the functions c*(A) it will follow that Ar < oo. Let 
Or = limAj,Ar c*(A), Oi G F, and Ay = Mr(ar). We assume that all Ay are distinct and 
define 

A^ = {Ar, rank(F) < R}. 

We assume that M-p has a finite number of critical points on [gmm, 5'max] for each F 
with rank(F) < R. Let , c^^ be all the local maxima of Mr. We assume that Mr(cf ) 
are distinct for all F with rank(F) < R and i. Define 

A2 = {Mr(cf), rank(F) < R, l<i< kr}. 

Let F be a cycle of rank r < R.T the cycle of rank r + 1 that contains F, and T a cycle 
that is contained in F\F. Let /r,T = {c : Mr(c) = Mt(c)}. We assume that the sets /r,T 
are finite and /ri,Ti H /r2,T2 = unless (Fi, Ti) = (F2, T2) or (Fi, Ti) = (T2, F2). Define 

A3 = {Mr(c), c G /r,T, rank(F) < i?, T C F \ F}. 

We assume that the numbers Mr (ax) are distinct for all choices of cycles F and T such 
that rank(F) < R, rank(T) < rank(F) and z/(F) G T. Define 

A^ = {Mr(ar), rank(F) < R, rank(T) < rank(F), i/(F) G T}. 

Finally, we assume that the sets A-^, A^, A^ and A*^ do not intersect and define 

A = {Ao, Ai, A^} := {0} U A^ U A^ U A^ U A^ U {00}, 

where we arrange A*; in the increasing order. 

Below we will define c*(A) on the successive intervals (Afe,Afc+i) using induction on 
k while assuming that A^ are known. The above definition of A^ and A*^ in terms of 
c*(A) does not constitute a circular argument, since we could instead define the pairs 
(Afe+i, c*(A) for A G (Afc, Afc+i)) inductively. Such an approach would lead to more compli- 
cated notations, though, so we avoid it. 
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Let us proceed with the inductive definition of c'(A). For A e (Aq, Ai) all cycles are 
passive and c*(A) = g{Oi) for all i. Assuming that the types of the cycles and the limits 
q{Oi) = limAiAfe c*(A) are known for A G {^k-i, Afc) with some < A; < m, we will describe 
the types of the cycles for A G (A^, A^+i) and specify the limits s(Oj) = limxix^c^{X). 
Then, assuming that the types of the cycles are specified for A G (A^, A^+i) and the values 
of s{Oi) are known, we will define the functions c*(A) for A G (A^, Afc+i). We distinguish 
a number of cases depending on whether A^ belongs to A^, A^, A^ or A^. 

First, however, we describe the procedure for determining the values of s(Oj) for Oj 
which belong to a cluster. 

Determining the values of s{Oi) and the types of cycles within a cluster. 
Suppose that we have defined s{Oi) — limAjA^ c*(A) for all that belong to a cycle F. 
Consider the cluster of cycles that are connected to F for A G (Afc_i, Afc). For each cycle 
F' in the cluster, we will define the values of s{0) for O G F' and specify its type for 
A G (Afc, Afc_|_i). 

First assume that i^(F') = Oj G F for A G (Afc_i, Afc). It will follow from the inductive 
construction that q{0') = q{0") if O', O" G F'. Let g(F') = q{0'). For O G F', we define 
s{0)=C{q{r),s{0,),X,,T'). 

For any cycle F" such that z^(F") G F', we can similarly determine the values of s{0) 
for O G F". Continuing this procedure inductively, we define the values of s{0) when 
O belongs to either of the cycles from the cluster. A cycle F' from the cluster will be 
engaged for A G (Afc, Afc+i) if Afc = Mt'{s{0)) for O G F' and active if Afc > Mt'{s{0)) for 

o G r. 

Case 1. Assume that Afc G Ai. Let F be such that Afc = A-p. For Oi G F, we 
define s(Oj) = C(g(Oi), g(i/(F)), Afc, F). The cycle will be engaged for A G (Afc,Afc+i) if 
Afc = Mr(s(0)) for O G F and active if Afc > Mr(s(0)) for O G F. 

The types of cycles that belong the the cluster connected to F for A G (Afc_i, Afc), and 
the values of s{Oj) for the equilibrium points in those cycles are determined according 
to the procedure described above. For the remaining equilibrium points O, we define 
s{0) = q{0). The remaining cycles don't change type. 

Case 2. Assume that Afc G A2. Let c be the local maximum of a cycle F such that 
My{c) — Afc. If F was not engaged for A G (Afc_i, Afc) or if q{0) ^ c for some O G F, then 
we define s{0) = qiO) for all the equilibrium points, and all the cycles have the same 
type on (Afc, Afc+i) as on (Afc_i, Afc). 

If F was engaged and qiO) = c for O G F, then for Oi G F, we define siOi) = 
C{q{Oi),q{u{T)),Xk,T). The cycle will be engaged for A G (Afc, Afc+i) if Afc = Mr(s(0)) 
for O G F and active if Afc > Mr(s(0)) for O G F. 

The types of cycles that belong the the cluster connected to F for A G (Afc_i, Afc), and 
the values of s{Oj) for the equilibrium points in those cycles are determined according 
to the procedure described above. For the remaining equilibrium points O, we define 
s{0) = q{0). The remaining cycles don't change type. 

Case 3. Assume that Afc G A3. Let F be a cycle of rank r < i?, F the cycle of rank 
r + 1 that contains F, and T a cycle that is contained in F \ F. Suppose that c is such 
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that Mr(c) = Mt(c) and = Mr(c). 

We define s[0) = q{0) for all the equilibrium points. All the cycles, other than 
perhaps F and T, have the same type on (Afc, A^+i) as on (Afc_i,Ajt). To determine the 
type of cycles T and T on (A^, Xk+i), we examine several cases. 

(a) If q{0) = c for all O e r U T, r and T were engaged, T was connected to F by 
a chain that contained only active cycles (other than T itself) and F was connected to T 
by a chain that contained only active cycles (other than F itself), then F and T becomes 
active. 

(b) If q{0) = c for all O e F U T, F was connected to T by a chain that contained 
only active cycles (other than F itself), but T was not connected to F by a chain that 
contained only active cycles (other than T itself), and T was not passive, then F becomes 
active on (A^., A^+i) if it was engaged on (Ajt_i, Xk) and becomes engaged if it was active. 

The type of T stays the same. 

(c) the same as (b) with F and T interchanged. 

(d) If none of the cases (a)-(c) applies, then F and T have the same types on (A^, A^+i) 
as on (Afc_i, Afc). 

Case 4. Assume that A^ G A4. Let A^ = Mr(aT), where cycles F and T are such 
that rank(F) < R, rank(T) < rank(F) and z^(F) G T. We define s{0) = qiO) for all the 
equilibrium points. All the cycles, other than perhaps F, have the same type on (A^, A^+i) 
as on (Afc_i, A^). 

The cycle F becomes active if it was engaged on (Afe_i, A^), q{0) — ar for all O G F 
and Mr{ar) < Ax- Otherwise, F has the same type on (Ajt, Xk+i) as on (Afc_i, A^). 

Now let us define the functions c*(A) on (A^, A^+i) assuming that the values of s{Oi) 
and the cycle types are known. For an equilibrium point Oj, we identify the cycle F with 
the smallest possible rank r such that Oj G F and the values of s{Oj), Oj G F, are not all 
the same. If no such cycle exists, that is if s{Oj), 1 < j < n, does not depend on j, then 
we define c*(A) = s(Oj) for A > Afc. 

Assuming that such a cycle F exists, let Fi, ...,F; be the cycles of rank r — 1 which 
comprise F, and let O G Fi. Here we number the cycles in such a way that A/'(Fi) = 
F2,...,A/'(F;) = Fi. Take the least j such that F^ is either passive or engaged (it can 
not happen that all the cycles Fi,...,Fi are active, since then all the values of s{0), 
O G F, would be the same, as follows from the inductive construction above). If F^ is 
passive, we define c*(A) = s{Oi) for A G {Xk,Xk+i). If F^ is engaged, we define c*(A) = 
C(r(Oi), C, A, Fj) for A G {Xk,Xk+i), where ( = +00 if Mr^. is locally increasing at r(Oj) 
and C = —00 if Mr- is locally decreasing at r{Oi). 

We can now summarize the above discussion. 

Theorem 5.1. Suppose that Assumption A holds and the hierarchy of cycles and the 
equilibrium points z/(F) for each cycle F of rank less than R do not depend on the choice 
of constants ci G [(71111117 fl'max] in the function a = ci{x , J^^^i CiXDi{x)) ■ Also suppose that 
the above assumptions on the sets A^, A^, A^ and A* hold. 
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Then the limits 

are uniform in x E Df for each 5 > 0, A G \ A, where the functions c*(A) were defined 
via the inductive procedure above. 



6 Example of a change in the hierarchy of cycles 

As in Section m we assume that there are three equihbrium points Oi, 02,03. For each 
ci,C2,C3 G [s'min, 5'max], the function /ci,c2,c3 is defined by ( fT9l) . We will assume that 
the hierarchy of cycles for a = a{x, fci,c2,c3{x)) depends only on C2. This is the case, for 
example, if d = 1 and Oi < O2 < O3. More precisely, suppose that there is c G ((y'min, 5'max) 
such that Assumption A holds for each choice of the constants q G [gmin, fl'max] such that 
C2 7^ c. We assume that Oi and O2 form a cycle T' = {Oi, O2} of rank one when C2 < c, 
while O2 and O3 form a cycle T" = {02,03} of rank one when C2 > c. 

As before, we will identify a number of "special" points and describe the asymptotic 
behavior of u'^(T^{\),x) for A G (A^, A^+i) and x G Df, i = 1,2, 3. In the process, we will 
make various assumptions about the quasi-potential that will be specific to the example 
at hand. 

In our example we assume that g{Oi) < (7(02) < c < (7(03). Define 



Ml2ic) = V;^[;;^l, M2l{c) = V;;^ Mp3(c) = V^,%'l, C G [grain, 



c 



M32{c) = M,,{c) = Mrn(c) = V^ir^\, M,r"{c) = V^%^^, c G (c,^7max]. 

Let Ai = Mi2{g{0i)) and A2 = M2i{g{02))- Define functions and by f |T3|) and 
(fT^ . respectively. Let A3 = inf{A : c^(A) > c^(A)}. Assume that at least one of the 
functions and is continuous at A3. Let c* = c^(A3) if is continuous at A3 and 
c* = c^(A3) otherwise. Let c^(A) = min(c^(A), c*) and c^(A) = max(c^(A), c*), A < A3. Let 
A4 = Mr/3(c*), As = sup,g[,.,^) Mp3(c) and Ag = M32(^(03)). 

Let us assume that Ai < A2 < A3 < A4 < A5 < Xq (see Figure 4). Define 



ci(A) = c2(A) = c^'(A) 



A3 < A < A4, 

min{c : c G [c*, c), Mr'3(c) = A}, A4 < A < A5, 



In order to formulate the results on the asymptotics of m'^(T'^(A), x) for A > A5, we need 
the functions d'^{X) and c^(A) defined as follows: 



d\X) = min{^(03), min{c : c G [c, giO^)], M^^ic) = A}}, A > A 



5; 



c3(A) = <! < A < Ae, 

' max{c, max{c : c G [c, (7(03)], M32(c) = A}}, A > Ag. 
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Figure 4: A case of three equilibrium when the hierarchy of cycles changes 



Let Ay = inf{A : d'^{\) > c'^(A)}. Assume that Ae < Ay and at least one of the functions 
d"^ and is continuous at Ay. Let c** = c?^(Ay) if d"^ is continuous at Ay and c** = c^(Ay) 
otherwise. Let As = Mr"i{c**) and assume that Ay < Ag. Define c^(A) = min((i^(A), c**), 
As < A < Ag, and c^(A) = max(c^(A), c**), < A < Ag. 
Let 

d\\) = min{c**,min{c : c G [c, c**], Mir"(c) = A}}, A > Ag, 

(A) = max{c, max{c : c G \c,c**], Mr"i{c) = A}}, A > Ag. 

Let Ag = inf{A : d^{\) > c^"(A)}. Assume that Ag < Ag and at least one of the functions d^ 
and is continuous at Ag. Let c*** = d^{Xg) if d^ is continuous at Ag and c*** = (Ag) 
otherwise. Define c^(A) = min(d^(A), c***), A5 < A and c^(A) = c^(A) = max(c'""(A), c***), 
As < A. 

Having thus defined the functions c*(A), i = 1,2,3, for all A > 0, we can now state 
that for each A > such that c* is continuous at A and every S > 0, the limit 

limM"(T"(A),x) = c\X) 

is uniform in x G -Df . 

On Figure 4, the limits lim^^o ^^(^^('^)5 2;), as functions of A, for x G Df, D2 and 
are depicted by thick, dotted and dashed lines, respectively. 
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